The quark core of protoneutron stars in the phase diagram of quark matter 
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We study the effect of neutrino trapping in new-born quark stars within a three-flavor Nambu- 
Jona-Lasinio (NJL) model with self-consistently calculated quark masses. The phase diagrams 
and equations of state for charge neutral quark matter in /3-equilibrium are presented, with and 
without trapped neutrinos. The compact star sequences for different neutrino untrapping scenarios 
are investigated and the energy release due to neutrino untrapping is found to be of the order of 
10 53 erg. We find that hot quark stars characterized, e.g., by an entropy per baryon of 1 — 2 and a 
lepton fraction of 0.4, as models for the cores of newborn protoneutron stars, are in the two-flavor 
color superconducting (2SC) state. High temperatures and/or neutrino chemical potentials disfavor 
configurations with a color-flavor-locked (CFL) phase. Stable quark star solutions with CFL cores 
exist only at low temperatures and neutrino chemical potentials. 

PACS numbers: 12.38.Aw, 12.39.-x, 24.85. +p, 26.60.+C, 97.60.-s 



I. INTRODUCTION 

The engine of explosive phenomena in astrophysics, 
such as gamma-ray bursts (GRBs) and type-II super- 
novae is presently not fully understood. The generation 
and propagation of neutrino fluxes, as well as the neu- 
trino interactions in the hot and dense nuclear matter 
envelope play key roles in models of the evolution of 
core-collapse supernova? (J). Detailed hydrodynamic sim- 
ulations of the gravitational collapse of a massive star 
show that the possibility to obtain a successful explo- 
sion depends to a large extent on the properties of the 
protoneutron star (PNS) that forms in the compressed 
core. About 99% of the gravitational binding energy is 
released by neutrino emission [1- 0- 13|- The supernova 
collapse and prompt neutrino production proceeds within 
milliseconds and the shock-compressed matter is heated 
up to about 30 — 50 MeV. At such high temperature the 
neutrino mean free path is much shorter than the radius 
of a PNS, R ~ 10 km, and neutrinos diffuse on a time 
scale of ~ 10 seconds. 

During this 'neutrino trapping' regime the number of 
neutrinos is quasi conserved and the neutrino chemical 
potential, fj, v , is of the order 200 MeV 0,0]. This state 
lasts until the temperature is low enough for the neutrino 
mean free path to become comparable to R. The behav- 
ior of the neutrino mean free path, which determines the 
timescale for the untrapping transition and the onset of 
the 'second shock' of the supernova, depends in a sensi- 
tive way on the microphysics of the hot and dense PNS 
interior. In particular, if superfluid and superconducting 
phases of hadronic and/or quark matter are created, the 
transport properties are affected and the consequences 
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for simulations of core-collapse supernova? could be im- 
portant. 

In a scenario where color superconducting quark mat- 
ter is the preferred state of matter at high density there 
are several new aspects to consider, for example: 

1. The kinetics of the phase transition, which even- 
tually proceeds with the release of latent heat and 
requires a nucleation time scale. 

2. The possibility that color superconducting phases 
have large pairing gaps (~ 100 MeV), which affect 
the transport properties in the PNS core. 

3. The possibility that color superconducting phases 
have a high critical temperature (oc the gaps) and 
therefore affect the formation and evolution of 
PNS. 

For recent reviews on dense color superconducting quark 
matter, see [l|,@,0,[l,[I|[l(J[ll|,[l2] and references therein. 

The energy released due to a phase transition to decon- 
fined quark matter in a PNS core can reach the order of 
100 bethe (=10 53 erg) 0, [3, [I3| , which is the correct or- 
der of magnitude for the energy of GRBs. Moreover, the 
nucleation timescale for a quark matter phase transition 
could explain the time delay statistics of GRB subpulse 
structure [l(| [jj}- It has also been emphasized that in 
the presence of a strong magnetic field, neutrinos propa- 
gating in hot superconducting quark matter can become 
collimated (beaming) and asymmetric, thus explaining a 
resulting kick velocity for the PNS [la ]. 

In the present work we consider a microscopic, al- 
beit schematic chiral quark model of NJL type, where 
the quark masses and pairing gaps are calculated self- 
consistently at the mean-field level, see [H, H3, Hj], H3, 
[23( 1 . We describe the phase structure of color super- 
conducting three-flavor quark matter for two different 
strengths of the phenomenological diquark pairing inter- 
action, in systems with and without trapped neutrinos. 
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The corresponding quark star solutions are considered as 
models for PNS cores and their properties are described. 
In particular, we are interested in the answers to the fol- 
lowing questions: 

1. What influence has the neutrino chemical potential 
on the phase diagram of quark matter? 

2. In which regions of the quark matter phase diagram 
can stable quark stars be found? 

3. Can both two-flavor superconductivity (2SC) and 
three-flavor color-flavor-locking (CFL) phases be 
realized under PNS conditions? 

4. How much energy can be released in the cooling 
and untrapping evolution for these models of a PNS 
core? 

The effect of neutrino trapping on the phase structure 
of color superconducting three-flavor quark matter has 
been investigated before in Refs. 0, HH, for homoge- 
neous phases, and in Ref. jH] for inhomogeneous phases 
of the LOFF type [27]. In this paper we focus on the 
stability of quark stars for different phase structures and 
estimates of the energy release in the cooling and untrap- 
ping transitions. 

The paper is organized as follows. In Section |Ul we 
define the NJL-type model that we use to describe color 
superconducting quark matter with trapped neutrinos. 
Section ITTT1 gives a summary of all results, which are dis- 
cussed in detail in subsections according to the main 
questions posed in the Introduction. The Conclusions 
summarize our main findings. 



II. MODEL OF HOT QUARK MATTER WITH 
TRAPPED NEUTRINOS 



Due to the high density in compact stars, strange 
quarks could exist in their interior. We therefore consider 
a grand canonical ensemble with up, down, and strange 
quark degrees of freedom. Because strange quarks have a 
relatively high mass, they should arise only at high den- 
sity. Therefore, at not too high densities, where matter is 
in a two-flavor state, there should be an excess of down 
quarks for matter to be charge neutral. Consequently, 
the difference between the chemical potentials of the up 
and down quarks can be sufficiently large for muons to 
be created by weak interactions, d <-» u + fi~ + v^. It 
is therefore necessary to include both electron and muon 
lepton degrees of freedom in the model. The r lepton 
is, however, not included, because it is too massive to 
play a significant role in compact stars. We neglect the 
influence of neutrino oscillations and therefore omit the 
r neutrino as well. 

The thermodynamics of the quark matter phase is de- 
scribed with an NJL-type model. The path-integral rep- 
resentation of the quark partition function is [20l . l2ll . [22j 



Z(T,A) 



VqDq e $ d ' lx ( c f+ c 9i+ c m) ^ 
q(i@ - rh + fi-f°)q, 



(la) 
(lb) 



C qq = G s^2(qT a q) 2 , (lc) 

a=0 

£ qq = G D ^2 (qij5TA^ACq T )(q T iCj 5 T A \ A q), (Id) 
4=2,5,7 

where fi and rh are the diagonal chemical potential and 
current quark mass matrices in color and flavor space. 
For a — 0, r = (2/3) 1 / 2 !/, otherwise r a and A a are 
Gell-Mann matrices acting in flavor and color space, re- 
spectively. C = i7 2 7° is the charge conjugation operator 



and q 



The coupling constants, Gs and Gd, 



determine the coupling strengths in the qq and qq chan- 
nels, which represent current-current interactions in the 
color singlet scalar meson channel and the scalar color an- 
titriplet diquark channel. We follow the argument in [2l[ 
concerning the tU(l) symmetry breaking in the pseu- 
doscalar isoscalar meson sector, which essentially is that 
the symmetry breaking is dominated by quantum fluctu- 
ations, and we therefore omit the 't Hooft determinant 
interaction. In the following we use the relative coupling 
strength 



V = G D /Gs 



(2) 



to parametrize the coupling in the diquark channel. In 
the choice of the four-fermion interaction channels we 
have omitted the pseudoscalar qq terms, which should 
be present in a chirally symmetric theory. These terms 
do not contribute to the thermodynamic properties of 
the deconfined quark matter phase at the mean-field 
(Hartree) level [2.3], to which we restrict the discussion 
in the present paper. 

After bosonization using Hubbard-Stratonovich trans- 
formations, we obtain an exact transformation of the 
original partition function ([I]). The transformed expres- 
sion constitutes the starting point for approximations, 
defined as truncations of the Taylor expanded action 
functional to different orders in the collective boson fields. 
In the following, we use the mean-field (MF) approxi- 
mation. This means that the bosonic functional inte- 
grals are omitted and the collective fields are fixed at the 
extremum of the action. The corresponding mean-field 
thermodynamic potential, from which all thermodynamic 
quantities can be derived, is given by 

^MF (T, p., flL e i M-Ljj ) 

1 ■ In Z MF (T, p) + fli (T, , fi L/t ) 
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Here, Mi — rrii + fa are the renormalized quark masses, 
rrii are the current quark masses, E a (p) are the eighteen 
independent quark quasiparticle dispersion relations, and 
fa (Ay) are chiral (diquark) gaps, see [2l[ for details. 
The gaps, fa and Ay, emerge from the auxiliary boson 
fields introduced by the Hubbard-Stratonovich transfor- 
mations and represent collective modes generated by qq 
and qq, respectively. f2j(T, , Rtu) i s the thermody- 
namic potential for an ideal gas of neutrinos, electrons 
(m e ~ 0.511 MeV), and muons (m M ~ 105.66 MeV) 
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where Ei = \Jp 2 + mf . fio is a divergent term that is 
subtracted in order to get zero pressure and energy den- 
sity in vacuum (i.e., at T = /i = 0) 



n = 



E 



- 6 



d 3 p 

W) 



^p 2 + {m l + $y. (5) 



The model has six conserved charge densities and asso- 
ciated chemical potentials. The U(3)f x SU(3) C x U(1)q 
symmetry of the quarks is broken in the presence of di- 
quarks. Therefore, there are only four mutually commut- 
ing conserved charge densities, e.g., the quark number 
density 



n = {q^q) = - 
two color charge densities 

n 3 = (q 1 hq) = - 
n s = (q 1 X 8 q) = - 
and the electric charge density 
riQ = (q^Qq) -n e - 



MF 



dfi 



MF 



dm 
SSImf 



dfl 



MF 



(6) 

(7a) 
(7b) 

(8) 



Here, Q — diag* (§, — 3, — 3) and n e (n M ) is the number 
density of electrons (muons). Consequently, the quark 
chemical potential matrix, fi, is 



fi = fj,+ Qn Q + A3/Z3 + A 8 /x 8 , 



(9) 



where fi is the quark number chemical potential, [iq the 
(positive) electric charge chemical potential, and /x 3 and 
Ms are color charge chemical potentials. Note the discus- 
sion of color neutrality in Ref. [2^ |. 



The quark number density is related to the baryon 
number density by 

n B =n/3. (10) 

The remaining two charges are the number densities of 
the lepton families 



n Lu =nn+ n Uu = - 



d^L c ' 



(11a) 
(lib) 



which are conserved when the neutrinos are trapped in 
the system and oscillations are neglected. Electrons and 
muons have both electric charge and lepton number, 
while neutrinos have lepton number only. Therefore 

(12a) 
(12b) 
(12c) 
(12d) 



Me = Mi e _ MQ; 

Mm = M£„ - MQ, 

Mfe = MieJ 



The lepton fractions, which represent the relative number 
of leptons and baryons, are defined as 



Yr 



Yr. 



n B 
n B 



(13a) 
(13b) 



Bulk matter in compact stars should be charge neu- 
tral. In NJL models there are no gauge fields that neu- 
tralize the color charge dynamically, because the gluons 
are replaced by effective point-like qq and qq interactions. 
Color neutrality must therefore be enforced by solving for 
the charge chemical potentials, /U3 and such that the 
corresponding charge densities, n a — (q^X a q), are zero. 
In addition, matter in compact stars should be electri- 
cally neutral and in /3-equilibrium with respect to weak 
interactions. The chemical potentials, /j,q, [13, and ^ 
are therefore determined such that the charge densities 
0-® vanish 

90mF 9&MF dtlyiF „ 



= 0. 



Observe that the definition of the quark chemical po- 
tential ® implies that matter is in /3-equilibrium. This 
holds true also when neutrinos are trapped and ^l, > 0, 
because, e.g., Me + = {ji Lt - Mq) + (-MzJ = _ MQ- 
The gaps, fa and Ay, are order parameters that are de- 
termined by minimization of the mean-field thermody- 
namic potential ® 



cT^MF 
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0, 



0. 



(15a) 
(15b) 



Local minima of the thermodynamic potential define 
competing phases and the global minimum is the physical 
solution. 
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At the charge neutral global minimum of the thermo- 
dynamic potential the pressure, entropy density, number 
densities, and energy density are 



P(T,/i,fi Le ,fi LlA ) = -SImf, 

dT ' 



n a (T, [i, fiL s , Hl) = - 



(16a) 
(16b) 

(16c) 



e(T, /j,, fi Le , fi Li J = -P + Ts + y^M g n a . (16d) 

The sum in the expression for the energy density should 
account for all conserved charges in the system. However, 
because we require that uq = 71.3 = n% = 0, only three 
terms are included. Once the current quark masses, cou- 
pling constants, and momentum regularization method 
has been fixed, equations ©-I©, ©, and - ((TSJ) de- 
fine a self-consistent set of equations for the quark matter 
model. 



A. Numerical methods 

The quark quasiparticle dispersion relations, E„.(p), in 
([3]) are eigenvalues of Hermitian matrices, see |2lj for de- 
tails. LAPACK is used to calculate the eigenvalues and 
the corresponding eigenvectors of these matrices. For 
given values of all parameters, the mean-field thermo- 
dynamic potential ([3]) is calculated with a Gaussian in- 
tegration quadrature. Observe that the integrands de- 
pend only on the magnitude of the three-momentum, so 
d 3 p — > Airp^dp. 

The derivatives of the thermodynamic potential are ex- 
plicitly calculated using the eigenvectors of the Hermitian 
matrices. The derivatives are of the form 



MF 



dx 



Fi(x) 



d 3 p 
J2i) 



3 



(17) 



for some simple functions F% and F2 that depend on the 
choice of the variable x. The numerical problem therefore 
is to calculate the derivatives of the quasiparticle disper- 
sion relations, E a (p), i.e., the derivatives of the eigenval- 
ues of the Hermitian matrices. It is easy to show that 



dE a _ X^X a 



x a x a 



(18) 



where X a is the eigenvector associated to an eigenvalue 
E a of a Hermitian matrix H . The derivatives of the Her- 
mitian matrices are sparse and all derivatives of the ther- 
modynamic potential are therefore obtained practically 
at the cost of computing the eigenvectors. The integrands 
of all derivatives and the thermodynamic potential are 
calculated in parallel and the momentum integrals are 
evaluated simultaneously to reduce computational redun- 
dancy. The gap and charge neutrality equations (fT4|l - lfT5|) 



are solved with a modified Newton method in multidi- 
mensions, i.e., essentially a steepest descent method, for 
different starting points in parameter space. The solution 
with the highest pressure is used. The results thereby ob- 
tained, and the consequences for the properties of quark 
stars are discussed in the next Section. 



III. RESULTS 

We use the same parametrization of the quark matter 
model as in 121] 



m u ,d 

771, 



G0A 2 



A = 



5.5 MeV, 
112.0 MeV, 
2.319, 
602.3 MeV. 



(19a) 
(19b) 
(19c) 
(19d) 



The relative diquark coupling strength, 77, is considered 
as a free parameter of the model. Here we present results 
for rj = 0.75 (intermediate coupling), which is the vacuum 
result obtained by a Fierz transformation, and 77 = 1.0 
(strong coupling) , which is motivated by the phenomenol- 
ogy of compact stars and heavy ion collisions [HI- Quark 
matter phases are characterized by the order parameters, 
4>i and Aij. In particular, the following phases have been 
identified in the numerical investigation of the model 

2SC phase: A us = A ds = and A ud ^ 0, 
uSC phase: A ds = and A ud ^ 0, A us ^ 0, 
dSC phase: A us = and A ud f 0, A ds ^ 0, 
CFL phase: A ud ^ 0, A us ^ 0, and A ds ^ 0. 

In addition, gapless phases exist, which are characterized 
by the presence of one or more quasiparticle dispersion 
relations that have no forbidden energy band above the 
Fermi surface. Such excitations exist when the differ- 
ences between the Fermi momenta and / or the renormal- 
ized masses of the paired quarks are sufficiently large [3C| . 
We denote gapless phases with a leading "g", e.g., g2SC 
for the gapless 2SC phase. 

In the following we present results for two different 
neutrino untrapping scenarios. In one scenario the ini- 
tial state is a hot, T = 40 MeV, quark core with trapped 
neutrinos, fj, Ve = 200 MeV. The core cools rapidly by neu- 
trino emission and a cold, T ~ 1 MeV, configuration with 
trapped neutrinos, /j„ = 200 MeV, forms. The mass de- 
fect due to cooling is obtained by comparing the masses of 
these two states for configurations of equal baryon num- 
ber. At low temperature, T ~ 1 MeV, the mean-free 
path of neutrinos increases and becomes comparable to 
the size of the core/star. The neutrinos therefore es- 
cape and the final state is a cold, T ~ 1 MeV, star with 
fjL v = 0. The mass defect due to neutrino untrapping 
is obtained by comparing the masses of [i u = 200 MeV, 
T ~ 1 MeV configurations with the masses of p, v = 0, 
T ~ 1 MeV configurations of equal baryon number. The 



5 



"final" state will continue to cool for millions of years with 
practically no effect on its structure, because tempera- 
tures below 1 MeV have negligible effects on the quark 
matter equation of state. Phase diagrams for interme- 
diate and strong coupling are provided for matter with 
trapped and untrapped neutrinos. The second scenario 
has a more "conventional" initial state, characterized by a 
fixed lepton number, Yt, e = 0.4, and entropy per baryon, 
s/ub — I — 2. In this case the temperature and neutrino 
chemical potential varies over the radius of the star. The 
properties of the quark star solutions and the effect of 
neutrino untrapping are similar in these two scenarios. 
While the initial state of the second scenario is moti- 
vated by detailed core-collapse simulations, it is more 
complicated and to some extent it provides less trans- 
parent results. 



B. Order parameters at fixed Yl c and s/na 

Next we consider a system with fixed values of the 
lepton fraction, Y Lc (T,^ u ) = 0.4, and the entropy per 
baryon, s(T, iiv)/nB(T, ji v ) = 1, 2. Consequently, for a 
given value of the quark number chemical potential, \i, 
the temperature and electron neutrino chemical potential 
are determined such that Yl s = 0.4 and s/ub = 1 (or 
2). These two equations are solved in parallel with the 
gap and charge neutrality equations. In Fig. [5] the tem- 
perature, constituent quark masses, charge chemical po- 
tentials, neutrino chemical potential, gaps, and number 
densities are plotted vs. the quark number chemical po- 
tential for charge neutral quark matter in /3-equilibrium 
at intermediate coupling, r\ = 0.75. Fig. [6] shows the 
same relationships for strong coupling, r) = 1. 



A. Order parameters at fixed fi 

For the calculation of quark matter phase diagrams in 
the plane of temperature and quark (or baryon) num- 
ber chemical potential, the gap equations for the order 
parameters of the model (the masses, Mi, and diquark 
gaps, A i: , ) are solved self-consistently. The values of the 
order parameters characterize the different quark matter 
phases, as described above. In addition to solving the 
gap equations, the color and electric charge neutrality 
conditions are enforced by solving for the corresponding 
chemical potentials (/jq, /J3, /is), for given values of the 
quark number chemical potential, /i, and temperature, T. 
In Figs.QJH the quark masses (upper row), diquark gaps 
(second row), chemical potentials (third row) and densi- 
ties of baryon number, electrons and muons (bottom row) 
are plotted vs. the temperature for three different values 
of the quark number chemical potential (columns) . The 
four Figures represent solutions without (fi u = 0) and 
with {^u = 200 MeV) trapped neutrinos for intermediate 
(tj = 0.75) and strong (rj = 1.0) coupling, respectively. 
Here we use the more compact notation 

\i v = pby e = fi Le , (20) 

for the electron neutrino chemical potential, because the 
muon neutrino chemical potential is zero as neutrino os- 
cillations are neglected. Note the gapless constraints, 
A? , that are included in the second row of these Figures. 
If Aij < Afj , the corresponding quasiparticle has a gap- 
less dispersion relation, i.e., when this condition is met 
there is no forbidden energy band above the Fermi sur- 
face. Typically, the gapless conditions are fulfilled only 
near the critical temperature (60 — 70 MeV for interme- 
diate coupling and 100 — 110 MeV for strong coupling), 
where the gaps go to zero in a second order phase transi- 
tion to a non-superconducting state. Since these temper- 
atures are well above the maximum temperature relevant 
for PNS, we can neglect gapless phases in our discussion. 



C. Equations of state 

The equation of state (EoS), e(-P), can be calculated 
using (fTS)) . e.g., for different values of the temperature, 
T, and lepton number chemical potential, fj,L s = H*v, or 
for fixed values of lepton fraction, Yl c (T, [i L J, and the 
entropy per baryon, s(T, /a^J/n^T, depending on 
the application. In Figs. [38] the EoS for three different 
pairs of (/J„,T) values are plotted for intermediate, r\ 
0.7"). and strong, 77 = 1.0, coupling, respectively. The 
energy density is discontinuous at the 2SC-CFL phase 
transition due to the discontinuity of the strange quark 
mass, which is obtained from the self-consistent solution 
of the system of gap equations lfl5|) . The phase transition 
is first order. When increasing T and/or the pressure 
at the transition increases. Fig. [9] shows the EoS for a 
fixed lepton fraction, Yl c (T, /j„) = 0.4, and two values 
of entropy per baryon, s(T, /j ly )/ns(r, /^y) = 1 and 2, 
for intermediate and strong coupling. When increasing 
the entropy per baryon, the pressure at the transition 
increases, whereas an increase of the coupling strength 
lowers it. 



D. Quark star sequences 

Given an EoS described in the previous Section, the 
corresponding compact star sequence is calculated by 
solving the Tolman-Oppenheimer-Volkoff equations for 
a static spherically symmetric object 

dP(r) _ [e(r) + P{r)][m(r) + 47rr 3 P(r)] 

7 — r 7\ 7 vi ' (^-La) 

dr r[r — 2m(r)\ 

m(r) = 4tt [ e(r')r' 2 dr', (21b) 

for different values of the central pressure. Due to the 
thermal pressure, the high-temperature EsoS do not ex- 
tend to zero pressure. The surface of hot configurations 
is therefore defined to be the point where the (approxi- 
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FIG. 1: (Color online) The constituent quark masses, gaps, charge chemical potentials, and number densities of charge neutral 
quark matter in /3-equilibrium with untrapped neutrinos, yL v = 0, at intermediate coupling, r\ = 0.75. The three columns (from 
left to right) represent solutions for /i = 400, 500, and 550 MeV. Ay are thresholds for the existence of gapless excitations 
[30l |. i.e., gapless excitations exist if Ay < A?-. 



7 



. M 

M -450 MeV 



M.-150MeV 



> 

s 50- 



M -100 MeV 



I , 1 , 1 , 1 , 1 o 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

20 40 60 80 20 40 60 80 20 40 60 80 

T [MeV] T [MeV] T [MeV] 




T [MeV] T [MeV] T [MeV] 




20 40 60 80 20 40 60 80 20 40 60 80 

T [MeV] T [MeV] T [MeV] 




T [MeV | T [MeV] T [MeV] 



FIG. 2: (Color online) The constituent quark masses, gaps, charge chemical potentials, and number densities of charge neutral 
quark matter in /3-equilibrium with trapped neutrinos, \l v = 200 MeV, at intermediate coupling, r\ = 0.75. Line styles as in 

Fig.m 
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FIG. 3: (Color online) The constituent quark masses, gaps, charge chemical potentials, and number densities of charge neutral 
quark matter in /3-equilibrium with untrapped neutrinos, fj, v = 0, at strong coupling, r\ = 1.0. Line styles as in Fig. [l] 



9 



M, 

M -450 MeV 



M -150 McV 



M, 

M -100 MeV 




50 100 50 100 50 100 

T[MeVl T [MeV] T |MeV] 

PIG. 4: (Color online) The constituent quark masses, gaps, charge chemical potentials, and number densities of charge neutral 
quark matter in /3-equilibrium with trapped neutrinos, /2 V = 200 MeV, at strong coupling, r] = 1.0. Line styles as in Fig. [T] 
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FIG. 5: (Color online) The temperature, constituent quark masses, charge chemical potentials, neutrino chemical potential, 
gaps, and number densities of charge neutral quark matter in /3-equilibrium at intermediate coupling (77 = 0.75) for fixed lepton 
fraction Yl c (T, (j, v ) = 0.4, and two values of the entropy per baryon, s(T , [i v ) / ub{T , = 1, 2. 




FIG. 6: (Color online) The same as Fig. [5]but at strong coupling, rj = 1. 
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FIG. 7: (Color online) The equation of state for charge neu- 
tral quark matter in /3-equilibrium at intermediate coupling, 
r\ = 0.75, and fixed values of the temperature and electron 
neutrino chemical potential. The discontinuities appear at 
the 2SC-CFL phase transition due to the Maxwell construc- 
tion, see text. 
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FIG. 8: (Color online) The equation of state for charge neutral 
quark matter in /3-equilibrium at strong coupling, r\ = 1.0, 
and fixed values of the temperature and electron neutrino 
chemical potential. 

mate) chiral symmetry is broken by a first order phase 
transition into the %SB phase. 

In Fig. [10] the quark star sequences for intermediate 
coupling, r) = 0.75, and fixed ((X U ,T) values are plotted. 
The discontinuity at the maximum mass configuration is 
a consequence of the appearance of CFL matter in the 
center of the stars. Stars with a CFL core are marginally 
stable and exist if T and [i v are not too high, see Fig. [141 
in the following Subsection. Increasing the coupling to 
i] = 1.0 leads to an interesting situation, shown in Fig.fTTl 




'o 500 1000 1500 2000 2500 
P [MeV/fm 3 ] 



FIG. 9: (Color online) The equation of state for charge neu- 
tral quark matter in /3-equilibrium for fixed values of the lep- 
ton fraction, Yl c (T, ^i„) = 0.4, and the entropy per baryon, 
s(T, /j,v)/n B (T, fj, v ). 



There is a small interval of star masses at M ~ 1.3 M Q 
for which mass twins occur: pure 2SC quark stars have 
stable, high-density mass isomers with a CFL quark core 
and a smaller radius. Fig. fTTl suggests an interesting sce- 
nario: upon mass accretion, a 2SC star could undergo a 
phase transition to its more compact twin with a CFL 
core, whereby binding energy is released. However, the 
presence of a hadronic crust could render the twin con- 
figurations unstable and this scenario therefore remains 
to be investigated thoroughly. In Fig. [12] we show the 
star sequences corresponding to the EsoS in Fig. [9] with 
fixed entropy per baryon (1 or 2) and fixed lepton frac- 
tion Yl b = 0.4. While the increase of coupling strength 
increases the maximum mass and the radius of the stars, 
both values of the entropy per baryon produce rather 
similar star sequences. A finite entropy per baryon of 
1 — 2 correspond to rather high temperatures that do not 
allow for stable CFL cores. In Fig. [13] the radial depen- 
dence of the temperature and baryon number density in 
PNS is illustrated for fixed values of the lepton fraction 
and entropy per baryon. These configurations, which re- 
semble realistic PNS formed in the adiabatic compres- 
sion of the cores in massive stars, have approximately 
constant temperature and neutrino chemical potential in 
their interior. Consequently, at the qualitative level the 
initial state of a PNS can be modeled by a fixed initial 
temperature, T(r) — T, and neutrino chemical potential, 
(J, v (r) = This result justifies the alternative neutrino 
trapping/untrapping scenario considered in this paper. 



E. Phase diagrams 

Next we present for the first time the central temper- 
atures and quark number chemical potentials of stable 




FIG. 10: (Color online) Quark star sequences for intermediate 
coupling, r\ — 0.75, and fixed values of the temperature and 
electron neutrino chemical potential. The discontinuity at the 
maximum mass configuration is a consequence of the appear- 
ance of CFL matter in the center of the stars. Stars with a 
CFL core are marginally stable and exist if the temperature 
and neutrino chemical potential are not too high, see Fig. 1141 
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FIG. 11: (Color online) Quark star sequences for the strong 
coupling, r\ = 1.0, and fixed values of the temperature and 
electron neutrino chemical potential. Stable stars with a CFL 
core exist if the temperature and neutrino chemical potential 
are not too high, see Fig. [15] 



quark star solutions in the phase diagrams of quark mat- 
ter, with and without neutrino trapping. In Figs. [141 
[TBI the phase diagrams of charge neutral quark matter 
in /3-equilibrium are shown for intermediate, f] = 0.75, 
and strong, r\ = 1, coupling and for two different values 
of the electron neutrino chemical potential, [i v = and 
200 MeV. 



FIG. 12: (Color online) Quark star sequences for fixed val- 
ues of the lepton fraction, Yl c (T, ^i„) = 0.4, and the entropy 
per baryon, s(T, fj, u )/riB{T, fi v ). Stars with a CFL core are 
unstable (marginally unstable for 77 = 1). 




r [km] 

FIG. 13: (Color online) The temperature and baryon number 
density vs. the coordinate radius of four different quark star 
solutions, for fixed values of the lepton fraction, Yl c (T, fit,) = 
0.4, and the entropy per baryon, s(T, /ub(T, p„). The 
mass of these quark star solutions is 1.25 Mq (1.4 Mq) for 
intermediate (strong) coupling. 
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FIG. 14: (Color online) Phase diagrams of charge neutral quark matter in /3-equilibrium at intermediate coupling, r\ = 0.75, 
for fixed values of the electron neutrino chemical potential, fi v = (left-hand side) and = 200 MeV (right-hand side). 
First-order phase transition boundaries are indicated by bold solid lines, while bold dashed lines represent second-order phase 
boundaries. The dotted lines indicate gapless phase boundaries and the thin solid lines are level curves of constant entropy per 
baryon. Hatched regions represent stable compact star solutions, with central quark number chemical potential fx(r = 0) = fi 
and temperature T(r = 0) = T. The cross-hatched regions correspond to baryon number twins, i.e., for these values of fi(r = 0) 
and T(r — 0) there exist stable 2SC stars and 2SC-CFL stars with equal baryon number. 




FIG. 15: (Color online) Phase diagrams of charge neutral quark matter in /3-equilibrium at strong coupling, r] = 1.0, for fixed 
values of the electron neutrino chemical potential, /i„ = (left-hand side) and \i v = 200 MeV (right-hand side). Line styles as 
in Fig. [ll 
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FIG. 16: (Color online) Radial phase structure of final states 
with CFL cores at strong coupling, r\ = 1.0. The T = 0, 
iiu = compact star sequence in Fig. [TT] has a branch of 
stable stars with CFL cores, this Figure shows the radial phase 
structure of these stars. 




H(r = 0) [MeV] 

FIG. 17: (Color online) Cooling of PNS with trapped neu- 
trinos, = 200 MeV, at intermediate coupling, 77 = 0.75. 
The dashed lines indicate the central temperature and quark 
number chemical potential of configurations with fixed baryon 
number, N = 0.2, 0.3, 0.4, . . . , 1.2 N . Solid lines indicate 
phase boundaries and level curves of constant entropy per 
baryon, see Fig. 1141 



As an upper limit estimate for the initial core temper- 
atures of quark stars we show the lines corresponding to 
an entropy per baryon of 3, which in the 2SC phase leads 
to a maximum temperature of ~ 40 MeV (~ 50 MeV) 
for intermediate (strong) coupling. For sufficiently high 
neutrino chemical potentials, quark stars with a CFL 
core are unstable. In particular, by comparing the left- 
hand (/i„ = 0) and right-hand (/i^ = 200 MeV) panels of 
Fig.[l5]it is clear that stars with a CFL core are rendered 
unstable by the finite neutrino chemical potential. For a 
typical value of /i„ = 200 MeV, which should be realized 
in the evolution of a PNS core [!, 5, we find no stable 
CFL cores. Note that the discussion of star temperatures 
exceeding the neutrino opacity temperature, T c ~ 1 MeV, 
makes sense only during the PNS era, i.e., when neutri- 
nos are trapped and the neutrino chemical potential is 
finite. Consequently, the region of stable CFL core stars 
shown in the panel on the left-hand side of Fig. [15] may 
not be realized at all during stages of hot PNS (quark 
star) evolution. We will return to this issue in the next 
Subsection. 

For temperatures below T c , we find a range of baryon 
numbers Nb = 1.39 — 1.43(5) N©, for which quark stars 
with a CFL core and a 2SC shell are stable. The phase 
structure of these stars is illustrated in Fig. HU Gapless 
phases of superconducting quark matter exist for both 
coupling strengths, at temperatures well above those of 
relevance for the cores of PNS, see Figs. fT4lfT5l The gap- 
less phases can therefore be neglected in the present dis- 
cussion of PNS evolution (3ll |. 



F. Cooling evolution and untrapping transition 

In this Subsection we consider a scenario of quark star 
cooling where the baryon number of the stars is con- 
served, i.e., we assume that there are neither accretion 
nor mass loss. In Figs. [T7lfT8l we show the cooling evolu- 
tion of PNS configurations in the phase diagrams of quark 
matter with trapped neutrinos, [i v — 200 MeV, for inter- 
mediate and strong coupling. The dashed lines indicate 
the central temperature and quark number chemical po- 
tential of configurations with fixed baryon number. The 
vertical solid lines indicate phase boundaries of the 2SC 
phase, while the horizontal lines denote curves of con- 
stant entropy per baryon. We observe that, as long as 
the neutrino chemical potential remains fixed, the cooling 
trajectories that start from stable configurations remain 
inside the region of stability. The reverse, however, does 
not hold. Upon heating, a configuration close to the CFL 
phase border may become unstable at high temperature. 
See, e.g., the trajectory for N = 1.6 N Q in Fig. [181 

Another instability occurs during the cooling evolution 
when a star reaches the neutrino opacity temperature, 
T c ~ 1 MeV, and neutrinos start leaving the system. 
The microscopic processes behind this neutrino untrap- 
ping transition are elastic and inelastic neutrino-quark 
collisions, which lead to a neutrino mean free path that 
exceeds the radius of the star as the temperature de- 
creases to T c . This characterizes the transition from 
the neutrino-diffusion regime with finite fi v to the free- 
streaming regime with //„ = 0, which takes place within 
the transport timescale of a few milliseconds [18]. Since 
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FIG. 18: (Color online) Cooling of PNS with trapped neu- 
trinos, \l v — 200 MeV, at strong coupling, r\ = 1.0. The 
dashed lines indicate the central temperature and quark num- 
ber chemical potential of configurations with fixed baryon 
number, N = 0.2, 0.3, 0.4, ...,1.6 N©. Solid lines indi- 
cate phase boundaries and level curves of constant entropy 
per baryon, see Fig. [15] 



this timescale is well above the strong and weak inter- 
action timescales, which establish the local thermal and 
chemical equilibrium, the untrapping transition can be 
considered as a quasi-static process that connects the 
initial state with finite neutrino chemical potential and 
the final state with vanishing fi v by a sequence of equi- 
librium star configurations. This justifies the estimate 
of the energy release due to neutrino untrapping from 
the mass defect between initial and final configurations. 
Cooling trajectories that are close to the CFL phase bor- 
der above T c continue within the CFL-twin region below 
it, i.e., conservation of baryon number requires that the 
phase structure of these stars change near T c . Leaving 
aside the details of the untrapping transition itself, we 
illustrate this situation in Figs. IT9ll20l for strong and in- 
termediate coupling, respectively. In these figures, the 
phase diagram above (below) T c represents quark matter 
with trapped (untrapped) neutrinos. We now discuss the 
strong coupling case illustrated in Fig. [19] in more detail. 
At T c the mean free path of the neutrinos is similar to 
the size of the stars and fi v — > 0. The central quark 
(or baryon) number chemical potential of configurations 
with fixed baryon number increases in the untrapping 
transition. When the neutrinos have escaped there are 
stable 2SC stars for N < 1.41 N and stable stars with a 
CFL core for 1.39 N < N < 1.43 N . The dashed lines 
in the CFL phase indicate the central temperature and 
quark number chemical potential of configurations with 
fixed baryon number, N — 1.4, 1.41, 1.42 N . The con- 
figurations with N = 1.4 and 1.41 N are CFL baryon 
number twins, i.e., for these baryon numbers there exist 



FIG. 19: (Color online) Neutrino untrapping transition at 
strong coupling, r\ = 1.0. The dashed lines in the 2SC phase 
(See Fig. I15|) indicate the central temperature and quark num- 
ber chemical potential of configurations with fixed baryon 
number, N = 0.2, 0.3, 0.4, . . . , 1.4 N . Above the critical 
temperature, T c ~ 1 MeV, of the neutrino untrapping transi- 
tion the neutrino mean free path is shorter than the size of the 
stars and [i v is therefore fixed at 200 MeV. This corresponds 
to cold PNS with trapped neutrinos. At T c the mean free 
path of the neutrinos is similar to the size of the stars and 
Hv — » 0. The central baryon number chemical potential of 
configurations with fixed baryon number increases in the un- 
trapping transition. When the neutrinos have escaped there 
are stable 2SC stars for N < 1.41 N and stable stars with 
a CFL core for 1.39 N < N < 1.43 N . The dashed lines 
in the CFL phase indicate the central temperature and quark 
number chemical potential of configurations with fixed baryon 
number, N = 1.4, 1.41, 1.42 N . The configurations with 
N = 1.4 and 1.41 N Q are CFL baryon number twins, i.e., for 
these baryon numbers there exist also pure 2SC stars. See 
Fig. US] and Fig. [2l]for further information. 



also pure 2SC stars. 

For intermediate coupling, shown in Fig. [20] there is 
only a narrow interval of baryon numbers for which a 
transition to marginally stable CFL baryon number twins 
is possible during untrapping. These baryon numbers are 
not shown in that Figure. 



G. Mass defects due to cooling and neutrino 
untrapping 

The cooling evolution of a hot PNS from T ~ 40 MeV 
to T c ~ 1 MeV and the subsequent neutrino untrapping 
transition entail significant changes to the EoS and the 
structure of the star. This results in a decrease of the 
gravitational mass. The mass defect, i.e., the difference 
of the gravitational masses before and after the cooling 
and/or untrapping transition, corresponds to the energy 
that is released from the star, predominantly by neutrino 
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FIG. 20: (Color online) Neutrino untrapping transition at 
intermediate coupling, r\ = 0.75. The dashed lines in- 
dicate the central temperature and quark number chemi- 
cal potential of configurations with fixed baryon number, 
N = 0.2, 0.3, 0.4, . . . , 1.1 Nq. CFL baryon number twins 
are marginally stable and exist only for a narrow interval of 
baryon numbers. They are therefore omitted. For further 
information, see Fig. [14] and Fig. 1191 



emission. A fraction of the neutrinos might, however, be 
converted into photons that give rise to a gamma- or X- 
ray burst. The PNS cool rapidly by neutrino emission. A 
few seconds after birth the stars are cold on the nuclear 
scale (T ~ 1 MeV) and further cooling has no direct im- 
pact on their structure. However, the neutrino mean-free 
path is sensitive to the temperature. Below some critical 
temperature (~MeV) the neutrinos that are trapped in 
the PNS escape and this has an effect on the structure. 
In Figs. 12111221 the energy release is estimated from the 
mass defects due to cooling and neutrino untrapping for 
strong coupling and two different values of the baryon 
number. 

In Figs. [23ll24l the energy release due to cooling and 
neutrino untrapping is plotted vs. the initial mass of the 
PNS, for intermediate and strong coupling, respectively. 
Fig. [25] shows the net energy release when PNS cool from 
a given entropy per baryon of 1 or 2 to zero and the 
neutrinos are untrapped, Yj, e = 0.4 — > 0. 



IV. CONCLUSIONS 

We have studied the effect of finite neutrino chem- 
ical potentials on the phase diagram and the finite- 
temperature EoS of three-flavor quark matter with self- 
consistently determined quark masses and pairing gaps. 
We confirm the results of Riister et al. [25|| that the 
phase transition to strange quark matter, such as the 
CFL phase is shifted to higher densities when neutri- 
nos are trapped in the system, thus making it unlikely 
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FIG. 21: (Color online) The mass defect due to cooling and 
neutrino untrapping at strong coupling, rj = 1.0. The dash- 
dotted line represents the sequence of hot (T = 40 MeV) 
PNS with trapped neutrinos (/x„ = 200 MeV). The PNS cool 
rapidly by neutrino emission. A few seconds after birth the 
stars are cold on the nuclear scale (T ~ MeV) and further 
cooling has no direct impact on their structure. However, 
the neutrino mean-free path is sensitive to the temperature. 
Below some critical temperature (~MeV) the neutrinos that 
are trapped in the PNS escape and this has an effect on the 
structure. The dashed (solid) line represent the sequence of 
cold stars with trapped (untrapped) neutrinos, i.e., T = and 
iXu = 200 MeV (nv = 0). When a hot isothermal 1.40 M PNS 
with baryon number N = 1.43 N© cools to low temperature 
by neutrino emission, the mass decreases to 1.35 M . At low 
temperature the neutrinos escape from the star and the mass 
decreases to 1.32 Mq. The relative mass defects are 3.5% due 
to cooling, and 2.1% due to neutrino untrapping. 



that strange matter exist in PNS cores. We find that in 
the presence of trapped neutrinos, the isospin mismatch 
induced by the /^-equilibrium conditions is reduced and 
the 2SC phase becomes more favorable. In particular, 
the onset of the 2SC phase is shifted to lower densities in 
the presence of trapped neutrinos. This result seems to 
be robust, as it has been found also for inhomogeneous 
phases of the LOFF type (26j. The investigation of the 
neutrino trapping effect on the LOFF phases in the phase 
diagram with self-consistently determined strange quark 
masses is an interesting task for future work. 

A new result of this study is the systematic analysis 
of the regions in the quark matter phase diagram that 
are realised in stable quark star cores. This analysis is 
carried out for two different strengths of the coupling in 
the diquark channel and for two different values of the 
neutrino chemical potential, i.e., with and without neu- 
trino trapping. A remarkable finding is that for the ini- 
tial evolution of a PNS, when the neutrinos are trapped, 
all configurations with a CFL core are unstable. In ab- 
sence of trapped neutrinos there are small, isolated re- 
gions in the temperature-density plane where stable con- 
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FIG. 22: (Color online) The mass defect due to cooling and 
neutrino untrapping at strong coupling, r\ — 1.0. Line styles 
as in Fig. [21] When a hot isothermal PNS with baryon num- 
ber iV = 1.4 N0 cools to low temperature by neutrino emis- 
sion, the mass decreases from 1.38 Mq to 1.33 Mq. At low 
temperature the neutrinos escape from the star and the mass 
decreases to 1.30 Mq. The final state is either a homogenous 
2SC star or a more tightly bound 2SC star with a CFL core. 
The mass difference between these two final states is 0.01%. 
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FIG. 24: (Color online) The mass defects due to cooling (T: 
40 MeV -> MeV) and neutrino untrapping (fx v : 200 MeV -> 0) 
at strong coupling, r\ = 1.0. 
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FIG. 25: (Color online) The total mass defects due to cooling 
(s(T, /ii,) /ub(T, /i„): 1 or 2 — > 0) and neutrino untrapping 
(Ylb{T, Hv): 0.4 — » 0) at intermediate, r\ = 0.75, and strong 
coupling, r\ = 1.0. 



FIG. 23: (Color online) The mass defects due to cooling (T: 
40 MeV — > MeV) and neutrino untrapping 200 MeV — > 0) 
at intermediate coupling, r\ = 0.75. 



figurations with a CFL core exists. A fraction of these 
2SC-CFL stars have baryon number twins in the 2SC 
state. We suspect, however, that the configurations with 
a CFL core would be unstable when the influence of a 
hadronic crust is taken into account. See, e.g., Ref. 
for an analysis of such hybrid star configurations. All 



2SC quark star solutions considered in this work are sta- 
ble, while the phase border in-between the 2SC and CFL 
phases marks the endpoint of stability for PNS sequences 
in the plane of central temperature and central quark 
number chemical potential. 

The maximum temperature of PNS are estimated to 
about 50 MeV for an entropy per baryon of 3, which is 
an upper limit for the entropy generated in the adiabatic 
compression of the core during collapse of a massive star. 
For this value of the entropy per baryon, the temperature 
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in the CFL phase could reach about 65 MeV, but the 
corresponding PNS configurations are unstable. 

When hot PNS cool by neutrino emission and reaches 
the neutrino opacity temperature, T c ~ 1 MeV, i.e., the 
temperature where the neutrino mean free path becomes 
comparable to the size of the quark star core, the neu- 
trinos created by electron capture in the collapse are un- 
trapped. By calculating and comparing the compact star 
sequences for different neutrino untrapping scenarios we 
find that the energy release due to neutrino untrapping 
(and/or cooling) is of the order of 100 bethe. This pro- 
cess could therefore be important for the inner engine of 
supernovas and gamma-ray bursts. 

The discussion of superconducting quark matter 
phases in the present work concerns the stability of PNS 
configurations with a CFL core and a neutrino heating 
mechanism for the outer core of newborn PNS due to 



the fact that the opacity temperature in (superconduct- 
ing) quark matter exceeds that in nuclear matter. A 
detailed study of the neutrino untrapping transition in 
hybrid stars is beyond scope of the present work and will 
be given elsewhere. 
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